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Abstract. Planar supersymmetric quantum mechanical systems with separable spectral 
problem in curvilinear coordinates are analyzed in full generality. We explicitly construct 
, the supersymmetric extension of the Euler/Pauli Hamiltonian describing the motion of a light 

O-f particle in the field of two heavy fixed Coulombian centers. We shall also show how the SUSY 

Kepler/Coulomb problem arises in two different limits of this problem: either, the two centers 
collapse in one center - a problem separable in polar coordinates -, or, one of the two centers 
files to infinity - to meet the Coulomb problem separable in parabolic coordinates. 



1. Introduction 

The classification of quantum Hamiltonians for which the Scrodinger equation is separable 
in curvilinear coordinates is well known. Eisenhart in Reference [Ij completely achieved this 
task in three dimensional systems. More recently Andrianov and loffe worked this problem in 
supersymmetric quantum mechanical systems both in two and three dimensions [2J. Our aim 
in this paper is to describe how SUSY quantum mechanical systems separable in polar, elliptic, 
and parabolic coordinates can be constructed in a rather geometrical way. We shall follow 
previous work of our group on the interplay of integrability and supersymmetry in some planar 
supersymmetric systems, both classical and quantal, see [31 HI El E]. In particular, we shall 
address a SUSY model built from the Euler two-center problem which is separable in elliptic 
coordinates. As a bonus, we shall show that the planar SUSY Kepler /Coulomb problem built 
by loffe and collaborators in [7J and independently by Wipf el al in [8] arises in two different 
limits corresponding respectively to the separability of the Kepler /Coulomb problem in polar or 
parabolic coordinates. 



2. M = 2 Supersymmetric Quantum Mechanics 
2.1. J\f = 2 Classical Supersymmetric mechanics 

We start by describing M = 2 SUSY classical systems such that the configuration space is 
a two-dimensional Riemannian manifold {M?,g^ji_y). Let q^, /u = 1,2, be a system of local 

coordinates in and "da, a, fj. = 1,2, a set of Grassman variables. A super-point in the 

configuration super-space C = M^'^ is given by: {q^,'da) = (g^, Q'^, "!?];, i?2i ^2)1 whereas the 
following commutation/anti-commutation rules hold: 

q^K-Kq^ = , <^^ + « = . 



The super-symmetric action may be expressed in the form (Einstein summation convention will 
be used along the paper): 

f ( I i 1 1 dW dW d'^W 

S = jdt ct + -9,. KDtK + -^R,.,r,mm - + ^a^^'^^ 

Here: W{q^,q'^) : R, is the supcrpotential that characterizes the interactions compatible 

with supersymmetry in the system. The covariant derivative for the Grassman variables is 

A^^ = ^?^ + ri;xC 

and Tup and Rnupr] snoe the standard Christoffel symbols and curvature tensor associated to the 
metric g. Choosing a zweibein ea in order to decompose the metric tensor g^^'^, i.e., gf^^ = e^,e^ 
curved and flat Grassman variables can be related in a simple way: i?a = 

The action functional is invariant with respect to the SUSY transformations 

S^q^ = e^'l , Si^'l = ier , Si^i^ = is (^'^'^^ + i^'^p^i^^' 

( dW \ 

52q'' = £#2' , '^2^?^ = -ie [g'^'-Q^ + , ^2^^ = 

where e is a Grassman parameter, and the Noether's theorem gives the conserved super-charges: 

dW dW 

Qi = 9,.rr, + —r,, Q2 = gp.q'r,-—r, (1) 

As a pre-quantization step we pass to the Hamiltonian formalism. A point in the reduced phase 
super-space is given by: {q^,Pu, i^i, ^2)- -'-^ ^^^^ space a super- Poisson structure is defined in the 
form: 

such that {p^,,q''}p = S""^, {p,i.,Pu}p = {q^,q''}p = 0, {Ct?^}p = ig^^Sap. It is easy to verify 
that the supercharges Qi and Q2 satisfy the classical SUSY algebra: {Qi, Qi}p = {Q2-, Q2}p = 
2iH , {Qa,H^P = {Qii Q2}p = 0, where H denotes the classical SUSY Hamiltonian function of 
the system. Finally, complex Grassman variables can be defined as follows: 1?^ = ^{^i^Ti'&i), 
and the corresponding complex supercharges are written as: 

Q+ = -^(gi-ig2) , Q- = ^iQi + iQ2) 

2.2. Quantization process 

Canonical quantization in the coordinate representation is performed by the promotions: 
{(f'^Pv) — > {q^,Pu), and ('0'^^, ■0'*). g'^ and p'^ become the operators: 

whereas the Fermi operators ijj^^' and ip^^'^ are defined in terms of the zweibein and the "planar" 
Fermi operators "0" and ip"""^: ijj'^ = Ca V'", '0'*^ = ea'ip^^ ■ The quantum operators q^^ and pi^ 
satisfy the commutation relations: 

[eM = iMi , \p^,p,] = = [e,r] , = 1,2 



and 'tj)^ and ■0'*^ verify the anti-commutation relations of the form: 

We win choose a representation for the Fermi operators in terms of the generators of the Chfford 
algebra of R^: 

= (7^ + il') , = (7^ + il') , V'^^ = (7^ - > V'^^ = (7^ - il') 

where 7^^, \Ji = 1, ... ,4, are the 4x4 Hermitian Euclidean gamma matrices: 
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cr% i = 1, 2, 3, being the Pauli matrices. Thus and ■0'*^ are explicitly written as: 
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and the quantum supercharges become 4 x 4-matrices of differential operators: 
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and oji and 002 denote the "spin connection contribution" to the supercharges. For orthogonal 
metrics, they read: 

The quantum super-Hamiltonian: H = ^{Q,Q^ = \ {QQ'^ + Q^Q^i is the 4x4 differential 
operator: 
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with a block-diagonal structure inherited from the eigen-spaces of the Fermi number operator: 
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The vacuum state |0) of the SUSY system is defined by the expressions: V^'^lO) =0, /i = 1,2, 
and thus the Fock space is built from this state as follows: The creation operators acting on |0) 
bring the system into one-particle (fermionic) states: ^''^^lO) = The two-particle (bosonic) 

state is: V'^'''|li) = II1I2) = — ^^^^112) = — |l2li)- A general state in this space, = Txiffi-^i© -^2, 
can be written linear combination of the form: 
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The super symmetric space of states is the direct product of J- with the Hilbert space L^(M^): 
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Thus Fermi operators can be viewed as the generators of the Clifford algebra of 
the space of four-component Euclidean spinors generated by 
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The supercharges move states between the different Fermionic sectors of the space of states: 

Q Q 



3. Separable SUSY Quantum Mechanical systems 

We shall consider J\f = 2 supersymmetric systems with Liouville models as the bosonic sector. 
The key property enjoyed by this class of systems is that they are separable using two- 
dimensional elliptic, polar, parabolic or cartesian systems of coordinates. Thus the target space 
for these systems is the Riemannian manifold equipped with the metric induced from the 
change of coordinates to two-dimensional elliptic, polar and parabolic coordinates, and the 
Euclidean metric for cartesian ones. All these coordinates are orthogonal, i.e.: g^^ = , V/x / z/, 
and thus there exists a "natural" choice of the zweibein: 
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The super-Hamiltonian ([3]) is written in these cases as follows: In the bosonic sectors the 
Hamiltonian acts by means of the scalar ordinary Schrodinger operators: 
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In the Fermionic sector, however, the super-Hamiltonian reduces to a 2 x 2-matrix Schrodinger 
operator with diagonals components: 
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and odd diagonal components: 
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3.1. Systems separable in elliptic coordinates: SUSY QM 

We first analyze the case of separability in elliptic coordinates. The Euler version of these 
coordinates is defined as: 
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J"^ . n and r2 denote the distances: 



and the inverse change is determined by: u 
r\ = \/ {x\ — d)^ + x\ , r2 = 1/ [x\ + d)^ + from a generic point (x^, x^) to two fixed points: 
(d, 0) and (— d, 0), the foci of the change of coordinates. Elliptic coordinates are restricted to 
the ranges: d < u < 00, —d < v < d. Thus this change applies the cartesian plane B?, with 
coordinates into the elliptic strip with coordinates q^ = u, q^ = v. The metric 

tensor is: 
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The bosonic operators in elliptic coordinates are written as: 
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and the natural zweibein: 
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determines the Fermi operators: ■i/'"^ = e^V'^^, V'" = e"'^"^, '^^^ = 62^'^^, ?/''^ = e2ip^. Finally the 
supercharges are expressed as 
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and allow to construct the super-Hamiltonian: He = \{Qe-,Q^e\'- 



(6) 
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3.2. Systems separable in polar coordinates: SUSY QM 

Reproducing the results for the case of standard polar coordinates: = r cos ip, = r simp, 
and identifying: = r,q'^ = (p, the Bosonic operators read: 

Pr = -^^^ ' r = r, p^ = -ih-^ , ip = ip, [f,pr] = [ip,p^] = ih 

The natural zweibein: e^^ = 1, = 0, = 0, = determine the Fermi operators: 
^r-f _ e^'ip^^, 'ip'^ = e\tl;^, 'ipf^ = 62 "0^^ and = 62 V'^, and the supercharges are written in 
polar coordinates as follows: 
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The general super-Hamiltonian is: 
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3.3. Systems separable in parabolic coordinates: SUSY QM 

Finally, we present the results corresponding to parabolic coordinates: = ^i, = ^2- The 

change of coordinates is given by: = ^ (^^ "^i); — ^1^2, with — oc < < oo, and 
< ^2 < oo. The metric tensor is: = (722 = Ci + ?2 ' 5i2 = 521 = 0, with a natural zwcibcin: 
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Thus the supercharges in parabolic coordinates are: 
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And the super-Hamiltonian is: 
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4. The SUSY two-fixed Coulomb centers problem 

We present as an example the super symmetric problem of two Coulombian fixed centers, 
determined by the quantum Hamiltonian: 

2m ri r2 

where the two centers are located at the points (— c?, 0) and (d, 0). ai and a2 denote the center 
strengths, and we assume, without loss of generality: ai = a > a2 = 5a, with 5 G (0,1]. 



Quantization of the charge imphes that: a = e'^Zi, e being the charge of the electron, and 
5 = 1^ < 1, with Zi, G N*, Zi > Z2. 

In order to construct a super symmetric version of this problem, we choose the superpotential 
of the system among the solutions of the Poisson equation: 

—V^W = V (10) 
2m 

in a similar way to what has been done in references [7] and [U [9] for the Kepler /Coulomb 
problem. We will follow the notation introduced in [5]. Equation (jlOp is separable in elliptic 
coordinates. Considering the ansatz: W{u,v) = F{u) + G{v), the equation is written as: 
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that reduces to a ODE system solved by the general expressions: 
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We shall restrict ourselves to the the simplest case where the separation constant k and the 
integration constants are zero: k = Ci = C2 = C3 = C4 = 0, 
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and thus expressions ([5]) and ([6]) for supercharges in elliptic coordinates become: 
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that, in turn, lead to the super-Hamiltonian of the system: 
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In [6] we have found many bound state eigenfunctions of Heo profiting from the separability of 
the corresponding Schrodinger equation in elliptic coordinates. Also, the partner eigenfunctions 
in the discrete spectrum of Hei have been constructed through the action of . We will 
present here only the specific results for the zero modes of the system. The zero energy 
wave functions for the scalar Hamiltonians can easily be found as those annihilated by the 
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Solving the corresponding differential equations, the solutions are: 
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where /„ and Kn represent the modified Bessel functions. Henceforth "^^^Ujv) is a ground 
state. We can write the zero modes in Cartesian coordinates: ^'q'^^(x^, x2) = S'^^\u,v), 

where S denotes the matrix of the change of coordinates in the 

spinor space: 
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The Fermionic zero modes are annihilated by the two supercharges: 
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and the resolution of this system of differential equations gives us two independent solutions: 
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Only the second one is normalizable, thus taking Ai = we find the Fermionic ground state: 
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^.i. The SUSY Kepler/Coulomb problems Limit 

There exists two natural ways to obtain the Kepler/Coulomb problem from the problem of 
two centers: 1. The two centers collapse in one, and 2. One center flies to infinity. Taking 
into account that polar and parabolical coordinates can be regarded as limiting cases of elliptic 
coordinates, we expect that the SUSY Kepler/Coulomb problem either in the polar separable 
version (Limit 1) or, in the parabolic separable version (Limit 2) will arise in these limits 
of the elliptic system of coordinates. If this is the case, a path to the explanation of the 
superintegrability of the Kepler /Coulomb problem will be open. 

Following [7] and [8] we consider as superpotential and super-Hamiltonian for the 
Kepler /Coulomb problem the following expressions: 
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for the 2x2 Fermionic block. 



4.I.I. The two centers collapse into one center. We will use the trigonometric version of the 
elliptic coordinates in order to guarantee the preservation of the orientation: u = cosh(^ — /?), 
V = cosry. r] G [0, 2tt], and /3 < ^ < oo is an arbitrary constant. 

{ X"*^ = dcosh(^ — /3) cos ?? , = (isinh(,^ — /?) sinr/ 

The limit from elliptic to polar coordinates is obtained by taking d = 2e^ and the arbitrary 
constant /3 — )• — oo. This limit produces d — )• and thus the two centers will collapse into only 
one. 

x^ = lim dcosh(^ — /3) cos r/ = r cos (/3 , = lim dsinh(^ — /3) sinr/ = rsinc^ 



where = r. Applying this hmits to the supercharges and the supersymmetric Hamiltonian of 
the two-center Euler/Pauh problem we obtain in a direct way the corresponding expressions for 
the Kepler /Coulomb problem: 
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with a = a{l + 6), i.e., the electric charge is the sum of the corresponding charges of the two 
centers. 

We have shown in [6] that eigenfunctions in the discrete spectrum of the bosonic super- 
Hamiltonian of the two center problem go smoothly to all the normalizable eigenfunctions in 
the Kepler /Coulomb system when the two centers collapse. The zero modes, however, behave 
in a different way such that an interesting situation arises. Taking the limit polar limit in the 
bosonic zero modes and its norm we find: 
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Thus, we check that the well-known bosonic zero mode of the SUSY Kepler /Coulomb problem 
expressed in polar coordinates arises in this limit. 

Simili modo, the Fermionic zero energy ground state of the SUSY Euler/Pauli problem 
becomes in this limit d ^ 0: 
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but the corresponding norm N diverges: 
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There is no normalizable Fermionic zero modes in the SUSY Kepler/Coulomb system. 

4-. 1.2. One center "escapes" to infinity. Using re-scaled parabolic coordinates, the Kepler/Cou- 
lomb superpotential and super Hamiltonian for the Kepler /Coulomb system read: 
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The limit from elliptic to parabolic coordinates not only requires to take d — )• oo, but also a 
translation of the origin: = x — d, x'^ = y, and a re-definition of the elliptic coordinates (u, v) 
as follows: 
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In the trigonometric version, these changes are written as: 6 = =b\/dsinh^ , 6 = \/d|sinT/|. 
Now, in the d — ?■ oo limit two-center operators are converted into the Kepler /Coulomb ones in 
parabolic coordinates: 

Q-p = hm Qe , Qp = lim > H-p = hm He 
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Moreover, the superpotential in this parabolic limit becomes the Euler/Coulomb superpotential 
in parabolic coordinates: 



M/(6,6) = lim Witv) = -^(C? + 

a— ^-oo n 

The Bosonic zero energy ground state of the SUSY Euler/Pauli problem goes in this limit to: 









i.e., to the zero mode of the SUSY Kepler/Coulomb problem in parabolic coordinates with 
electric charge a. 

The Fermionic zero energy ground state of Euler/Coulomb problem goes to a zero energy 
eigenstate of the Kepler /Coulomb super Hamiltonian: 



V 






ma ( c2 t2\ 





but the norm of this state is again infinite. There are no normalizable Fermionic zero modes of 
the SUSY Kepler problem both in the polar and parabolic regimes. 
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